THE LIFSHITZ SLYOZOV WAGNER EQUATION 
FOR REACTION-CONTROLLED KINETICS 



APOSTOLOS DAMIALIS 



Abstract. We rigorously derive a weak form of the Lifshitz-Slyozov— Wagner 
equation as the homogenization limit of a Stefan-type problem describing 
reaction-controlled coarsening of a large number of small spherical particles. 
Moreover, we deduce that the effective mean-field description holds true in the 
particular limit of vanishing surface-area density of particles. 



1. Introduction 

The latc-stagc behavior of a material undergoing a first-order phase transition 
(due to changes in temperature and/or pressure for example) is characterized by 
thermodynamic instability resolved through phase separation and consequent coars- 
ening of the emerging phase. In the case of the new phase occupying much smaller 
volume fraction, and thus appearing as well-separated particles, this coarsening pro- 
cess (known as Ostwald ripening) is driven by the minimization of surface energy at 
the interface via diffusional mass exchange between particles while the total mass 
or volume of each phase is conserved. The result of this kind of mass diffusion from 
regions of high to regions of low interfacial curvature is the growth of large parti- 
cles and the shrinkage and final extinction of smaller ones. For a review of some 
aspects of Ostwald ripening, mainly from the physical and modeling viewpoint, see 
the survey by Voorhees [21] or the book by Ratke and Voorhees [18] . 

In this coarsening scenario the mass-diffusion process can be controlled by two 
different mechanisms: either by the diffusion of atoms away from the particles and 
into the bulk, or by the reaction-rate of attachment of atoms at the phase interface. 
In the former case (diffusion control), the random exchange of atoms between the 
particles and the bulk is sufficiently rapid and the surrounding of each particle is in 
thermal equilibrium with the atoms in it; in the latter (interface-reaction control), 
detachment and attachment are slow compared to diffusion and the surrounding 
bulk can be out of equilibrium with the particle interface. We refer to the physics 
literature for more details, for example, Slezov and Sagalovich [19) . Bartelt, Thcis, 
and Tromp [3]; for a related mathematical treatment see Dai and Pego [5]. 

The classical theory for Ostwald ripening was developed by Lifshitz and Slyozov 
[2] and Wagner [55] in the case of supersaturated solid solutions in three dimensions. 
The Lifshitz-Slyozov- Wagner theory statistically characterizes the evolution by the 
particle-radius density n(t,R), where n(t,R)dR is defined to be the number of 
particles with radii between R and di? at time t per unit volume. In the late 
stages of the phase transition nucleation and coalesence of particles can be neglected 
since new nuclei dissolve immediately and since particles cannot merge because of 
the large distances between them. Thus, the particle-radius density satisfies the 
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continuity equation (see [T8l §5.1]) 

-n(t,R) + —(v(t,R)n(t,R))=0, 

where v(t, R) denotes the growth rate of particles of radius R at time t. Using a 
mean- field ansatz (cf. Section [3]), Lifshitz, Slyozov, and Wagner formally calculate 
that 

W n(t, R) + — [^{Ru l)n(t, R)J = 0, 

with 

poo I poo 

u(t) = J n(t, R)dR J Rn(t, R) dR, 
in the diffusion-controlled case, and 



ith 



pOO I pOO 

u(t) = J Rn(t,R)dR J R 2 n(t,R)dR, 



(1.1) 



in the reaction-controlled one. both results valid in the limit of vanishing mass or 
volume fraction of particles. 

In [TT] and [T3] Niethammer rigorously derived the effective equations in the dif- 
fusion-controlled case, starting from a quasi-static one-phase Stefan problem with 
surface tension and kinetic undercooling, 

-Au = in Q \ G, 

V = Vu ■ n on dG, 

u = H + PV on dG, 

and restricting it to spherical particles. The same was also done in [TT] for the full 
time-dependent parabolic problem but without the kinetic-drag term (3V . Here, u 
is a chemical potential, n is the outer normal to the particle phase G, V is the nor- 
mal velocity of the phase interface dG, and H is its mean curvature. The domain 
SI C M 3 is considered bounded and (3 is a parameter that comes from the nondi- 
mensionalization and scales like diffusivity over mobility. The second boundary 
condition is the Gibbs-Thomson law, coupling the curvature of the interface with 
the chemical potential, modified by accounting for kinetic drag. Note that while 
under diffusion control the parameter f3 is small and the kinetic drag can even 
be neglected (thus yielding the well-known Mullins-Sckerka model [ID]), in the 
reaction-controlled case the values of j3 are large and, therefore, the kinetic-drag 
term is necessary. For a derivation of such sharp-interface free-boundary problems 
from continuum mechanics and thermodynamics see the book of Gurtin [S] . 

The goal in the following is to use the techniques developed in [TT] and [12] to 
derive the effective equations in the reaction-controlled case. This involves passing 
over to a different time scale incorporating the parameter (3 tending to infinity 
(see Section [2]) and, as a result, some extra manipulations in the proofs. Except 
for the scaling, in Section [5] we also give short proofs of some useful preliminaries 
and discuss the validity of the mean-field description while in Section [3] we prove 
pointwise estimates for approximate solutions and for the growth rates of particles. 
Finally, using these estimates, in Section [4] we pass to the homogenization limit of 
infinitely-many particles and obtain a weak form of the Lifshitz-Slyozov-Wagncr 
equation. 
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In comparison with the results in the diffusion-controlled case, we make precise 
that the crucial quantity that has to vanish in order to neglect direct interactions 
between particles and justify the expected mean-field law is the surfacc-arca density 
of the particles in contrast to their capacity in the other case (see [TT] and |12j). 
This difference is of interest since the asymptotic limits of vanishing surface area 
and capacity have different physical interpretations and further refine the naive 
general limit of vanishing mass or volume. For the reaction-controlled case though, 
the result is in some sense to be expected since the limit of vanishing surface- 
area density corresponds to the physics of the interface-reaction-controlled scenario, 
where there is an obvious dependence on the area of the interface. 

2. Formulation, scaling, and preliminary estimates 

We start with problem (| 1 . 1 [) where the quasi-static approximation to the para- 
bolic diffusion equation is justified by the small intcrfacial velocities present during 
late-stage coarsening. (See the discussion in Mullins and Sckcrka [TO].) 

We further suppose that the solid phase consists of spherical particles with cen- 
ters fixed in space, a simplification that can be justified by the work of Alikakos 
and Fusco pQ, [2], and Velazquez [20] . Denoting these particles as Bi, where each 
Bi is the closed ball B(xi, Ri(t)), the particle phase is then the union UBi and its 
isotropic evolution can be modeled by averaging the flux in the Stefan condition, 
i.e., 

V = Ri(t) ■■= i Vu-n, 

J dBi 

where the average integral is defined as 

£ f:= w\L f ' 

for a function / on some domain D, and where the overdot denotes a derivative 
with respect to time; the Gibbs-Thomson law becomes then 

u = — +0Ri, 

Hi 

since in the case of spheres the mean curvature is the inverse radius. 

To have many small particles in a bounded domain, for a system with size of order 
0(1), say the unit cube [0, l] 3 , let 8 be the typical particle distance with < 8 -C 1. 
For the distribution of particle centers in space, we assume, for simplicity, that they 
are situated on a three-dimensional lattice of spacing 8. Then, the initial number 
density of particles Ni(8) will be bounded by l/<5 3 , and for the particles to be small 
let the typical particle size be 8 a for a > 1. For times t £ [0, T] we choose a 8 small 
enough so that adjacent particles of size 8 a will not collide during the evolution up 
to a maximal time T. 

Concerning the assumption on the spatial distribution of particles, a more general 
assumption like inf,^j \xi —Xj\ > c8, for a constant c > 0, would still be enough for 
our purposes in this work. These considerations will also be used in the proof of 
Lemma [3 . 2 1 where we approximate a certain sum over all particles by an integral. For 
an approach using more sophisticated deterministic and stochastic assumptions on 
the distribution of particles with respect to homogenization we refer to Niethammer 
and Velazquez [15], [T7], where also further refinements of the theory are made. 

To have particle sizes of order 0(1) as well, we rescale 



4 



APOSTOLOS DAMIALIS 



and motivated by the scaling invariance of problem (cf . [5] ) , 

u s •= S a u t 5 ■= — 

Notice that this rescaling is another way of addressing the reaction-controlled 
regime. Instead of rescaling time by (3 and then letting f3 tend to infinity, we 
keep P fixed and positive, and specially rescale as above letting 6 tend to zero. 
Since now (3 plays no significant role, we will set it to unity in what follows. In 
addition, one easily sees that the transformations Rf , u s , and t s preserve the form 
of the equations. From hereon we also drop the superscript S from the notation for 
time and to denote the dependence on the new scale we write 



Bf := B( Xi ,6<*Rf). 

Finally, note that under diffusion control the relevant scale for time would be i5 3a 
instead of S 2a . This difference is key to all that follows, leading to different consid- 
erations on the validity of the mean-field model. (Cf. the remarks following Lemma 

O) 

As initial data, for every particle-center xi we associate a corresponding bounded 
initial radius Rf (0) with the assumption that 

su PieJVi i?f (0) < Ro, 

uniformly for some constant Rq- To consider a closed system, we impose a no-flux 
Neumann boundary condition on the outer boundary of SI, i.e., 

Vu s • n = ondil. 

In case the ith particle vanishes at time ti := sup{t | Rf (t) > 0}, for times later 
than ti we define Rf to be zero, reduce the number N(t) := {j \ Rj(t) > 0} of 
active particles by one, and neglect the boundary dBf in the boundary conditions. 
In the following, all sums, unions, and suprema will run over the set N(t), with 
iV(0) = iVi, and any further reference to the particle- number density will mean the 
active particle-number density TV unless otherwise noted. 

Summarizing, the restricted and rescalcd problem for the particle radii can be 
considered as a nonlocal, TV-dimensional system of ordinary differential equations 

i t _ , 



for times t £ (0,ti), ti < T, and with bounded initial data Rf{0) for every i, while 
the chemical potential is determined by 

-Au s (t,x) = in n\UjBf(t), (2.2) 
u 5 (t,x) = j^+Rf(t) ondBf(t), (2.3) 

and the Neumann condition on the outer boundary. 

Global existence and uniqueness of continuous, piecewise-smooth solutions for 
a similar restricted Stefan problem was proved in [12j by an application of the 
Picard-Lindelof theorem, the only difference being the different time scale. These 
solutions are not globally smooth due to the singularities arising from the extinction 
of particles; however, they are smooth in the intervals between the extinction times 
ti. In the following, when we mention solutions of the problem we will mean such 
continuous, piecewise-smooth solutions that exist up to any given time T. 
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It is easy to see that equations (|2.ip . (|2.2[) . (|2.3p . along with the outer boundary 
condition conserve the volume and decrease the interfacial area of the particle phase. 
Indeed, differentiating the total volume of particles with respect to time gives 

j- y RUt) 3 = 3 y Rt{tfM{t) = 3 y nf®* , n * „ / w-n 

where the last sum vanishes due to the divergence theorem, equation (|2.2p . and the 
no- flux condition on dtt. The decrease of total surface area follows from the next 
a priori estimate. 

Lemma 2.1. For any time t £ (0, T), the solutions of the problem satisfy the 
following energy equality. 

Proof. Multiplying —Au s = with u s , integrating over \ UBf, and integrating 
by parts gives 

/ \Vu s \ 2 + y f (Vu s -n)u s -[ {Wu s -n)u S = 0, 
Jn\uBf i JdBf Jan 

where the last term vanishes due to the Neumann condition on the outer boundary. 
Thus, using equations (|2.3p and (|2.ip we get 



i 

and after rearranging, 

]T(i?f) 2 |Af| 2 +£i?fAf + -^ / |V U 5 | 2 =0. (2.4) 

The result follows from an integration over time. □ 

After normalization with respect to the initial particle- number density Ni, this 
energy equality can yield useful information on the validity of the mean-field ap- 
proach. In fact, we have 



where the right-hand side is uniformly bounded by the assumption on the initial 
radii. For the left-hand side to stay bounded as well, if the quantity NiS 2a tends to 
zero, the same must hold for \S7u s \ and it is exactly this limit of vanishing surface- 
area density of particles that results in a mean field that is constant in space since, 
in particular, 

\7u s -> in L^TjiT^fi)). 
Here and in the following, to obtain global estimates that are uniform in S we extend 
u s to the interior of particles, and thus to the whole of il, by its boundary values. 
It is important to note that in our scaling setup, for the surface area to vanish as 5 
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tends to zero, the exponent a must be strictly larger than 3/2 since Ni is 0(l/<5 3 ). 
These facts will be made precise in Corollary 13.31 where we give an estimate of the 
mean- field effect. Note also that we do not address here the critical case a = 3/2 
that corresponds to finite surface area. For that one would have to use the different 
methods developed by Niethammer and Otto in |13j . 

Finally note that for similar considerations under diffusion control, the corre- 
sponding quantity would be the capacity NiS a due to the different time scale. In 
three dimensions, this capacity effect fits to general homogenization results as in 
the work of Cioranescu and Murat [4] ; to our knowledge though, the surface-area 
effect has not been explicitly discussed in the relevant literature. 

3. Approximation and growth- rate estimates 

As in the mean-field ansatz of Lifshitz, Slyozov, and Wagner, we suppose that 
the system is dilute enough so that particles behave as if they were isolated and we 
base our approximation on the solution of a single-particle problem. 

Consider problem (|2.1|) . (|2.2|) . (|2.3p for a single spherical particle centered at 
the origin and with initially unsealed radius r that we rescale as r d := r/5 a , along 
with the corresponding reaction-controlled rescalings for a chemical potential u 5 r 
and time, as in Section [2l For this rescaled particle Bf we consider the following 
problem in the whole space: 

**(*) - JZhL?^ I VuJ-n ondB s r , 



4,TrS 2a r s (t) 2 J dB s 
where the chemical potential uf(t,x) satisfies 

-Au*(t,x)=0 inR 3 \B*, 

u 5 r (t,x) = + r s (t) forxedB 5 r , 

and the mean-field assumption is posed as a condition at infinity, i.e. 

lim uf.(t,x) = u s r (t). 

\x\ — >oo 

This problem can be explicitly solved to give 

i« r $ (f\ Xa 



and 



1 + S a r s (t) V r s {t) t 
Note that in the formal limit of 5 tending to zero, the expected effective equations 
take the general form 

u(t,x)=u(t) and r~u , 

r 

as in the reaction-controlled Lifshitz-Slyozov- Wagner theory. 

Going now back to the many-particle problem, a calculation using the single- 
particle growth rate above along with the requirement that the volume is conserved 
gives the following expression for the mean field 

fi'=V /V ( fl '"> 2 (31) 
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The effect of this mean field plus a sum of single-particle solutions will be the 
monopole approximation to the solution u s supposing that there arc no direct 
interactions between particles. To this end, let us define the approximate solution 



( s (t,x) := u s (t) + Y 5 "?*®J l - u S (t)Rf(t))-^ 
sv ' ' ' 4^ l + S a R 5 t {t) K y ' lK> '\x-Xi\ 



(3.2) 



for x G n\UBf{t). 

Below is a maximum principle tailored to our setting that will be used to compare 
the approximation and the solution in the lemma next. Its proof can be found in 

Lemma 3.1. Let Q be a Lipschitz domain and let LiBj d il be a finite collection of 
disjoint closed balls. Then, a function v which is constant on each of the boundaries 
dBi and satisfies 

-Av = infl\UB t , 



Vw • n > on dBi 



I dBi 

Vv ■ n > on dil, 
where Cj > for all i, also satisfies 

v>0 infl\UBi. 

Lemma 3.2. For any time t G (0, T) and small positive s, the chemical potential 
and its approximation satisfy 



II" - C lli°°(n\uB*)(*) < c6 



2a-3-e 



sup Rf(t) (I + u\t) sup Rftt)). 



Proof. Since the difference u s — ( 5 is already harmonic in fi \ UBf as C 5 is a su- 
perposition of fundamental solutions, we would like to estimate to what extent it 
satisfies the maximum principle's boundary conditions. 

For the condition on the particle boundaries, we use equations (|2.1[) . (|2.3p . and 
the definition of ( s to calculate for x on the boundary dBf of the ith particle, 



C\t,x) 



VC • n - u 



%*)-{ 

Jdl 



VC • n 



dB" 



1 r 5 a R S i , . , 



5° 



' 3 \\X-Xj\ 

and since by the divergence theorem there holds for j ^ i, 

5° 



V 



5 a 



\X-Xj\ 



while for j 



{ 

f V 

JdBf 



\X-Xj\ 



n = 0, 



5"(Rf) 



S\2 ■ 
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we continue the calculation to get 
1 1 - u s R$ 



u s - 



l + 5 a RV 3 ' 

t L \ J ** 



5° 



\X-Xj\ 



S a R 



5° 



S a Ry 3 '\x-Xj\ 



< 5 2a ~ 3 sup R s j(l + u s sup Rj)J2 



6 3 



< C5 2a ~ 6 sup R] (1 + u 6 sup R°). 



(3.3) 



In the last step, keeping in mind the assumptions on the spatial distribution of 
particle centers, the sum is bounded for j ^ i since it is considered as a Riemann- 
sum approximation to the integral 

, — ^ — r dy, 

a\x-y\ 

which in turn is bounded using radial symmetry around the singularity and where 
the factor S 3 in the sum compensates for the scaling in space. 

To further fulfil the maximum principle's outer boundary condition on <9f2, we 
consider the comparison function ( s + z s , where the auxiliary function z s solves the 
problem 



VC • n in CI, 



IdQ. 

Vz s ■ n = -VC 5 • n on dVL, 



0. 



(3.4) 



such that the comparison function C 5 + z 5 has zero normal derivative on dil. To 
work with the maximum principle, z s also needs to be harmonic in ft and for that 
we need that the integral J gn VC 5 • n vanishes. But, 

/ VC* • n = 5 a y Sa *^ (1 - R*u 5 ) [ V(V^- 
Joa V! + M r Jen \\x-x 



n. 



where the last integral equals — An, independent of i. Thus, z is harmonic if and 
only if 

Rf (Rf) 2 



? l + 5 a R 5 J ? 1 



S a R: 



5 > 



which is exactly the mean field (|3.ip as dictated by the single-particle ansatz in 
the beginning of the section. Moreover, since now z s is harmonic, the divergence 
theorem further gives 

/ Wz s -n = 0. 

JdBf 

A construction as in Lemma 3 of and elliptic regularity theory (see Gilbarg 
and Trudinger [7]) give the estimate 

||**IU~(n) < C E S 2a - 3 - £ supRf(l + u s supRt) : 
where £ is a small positive number. 
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Let us now apply the maximum principle to the function 

/+ := u s - C S -z s + CS 2a - 3 - £ sup Rf (1 + u 5 sup ijf ). 

For a large enough constant C, say 2C e , the following hold for /+: it is harmonic, 
there holds V/+ • n = on dil by the construction of z 5 , and for the constants 
Cj = l/47T(5 2Q (i?f) 2 , estimate (jOJl gives 

- C 5 - z s + CS 2a - 3 - £ supiZf (1 + u 5 sup Rt) - c, [ S7(u s - ( s - z s ) ■ n > 0. 

Thus, / + satisfies the maximum principle's conditions and therefore, /+ > in 
n\UBf, i.e., 

u 5 - C 5 - z" 5 > -C5 2a - 3 - e svpRj (1 + u" 5 sup-Rf). 

Using the maximum principle with — u instead of v, the function 

/_ := u s -C s -z s - C6 2a - 3 - £ supRf (1 + supi?f), 

again satisfies the corresponding conditions and, as above, yields /_ < in Sl\UBf , 
i.e., 

u s - C 5 - z s < CS 2a ~ 3 - £ supflf (1 + S A supi?f). 
Combining the last two inequalities, we get 

||u 4 - C 6 ~ A\l°°(q\ub*) < C6 2a ~ 3 - £ sup Rf(l + u 5 su P Rf) 
and the lemma follows by the triangle inequality using the regularity of z s . □ 

In the previous lemma it is clear that our approach excludes the critical case 
a = 3/2. In the following we introduce, for technical reasons, a new exponent 
7 > with the property 

<F := m&x{5 a ,6 2a - 3 ,5 2a - 3 - £ } 

for each a greater than 3/2 + e. 

As a corollary to the previous lemma we can now estimate the effect of the mean 
field. 

Corollary 3.3. For any time t £ (0,T) and 7 > 0, the chemical potential and the 
mean field satisfy 

\\u s - u fi || L - ( n\uBf)(*) < C<5 7 (1 + 2supitf(f))(l + u 5 (t) sup Rf(t)). 

Proof. By the triangle inequality and Lemma 13.21 there holds 



W S -^\\L^<n\uB^)<\ I C 5 ~ u 5 \\l~ rrnuB?) + CS 2a ~ 3 - £ sup B* ( 1 + u 5 sup R s 



L°°(fi\U.B?) ^ IIS - " HL°°(f2\UB?) ~r ^ u P ^jv ' aayrLjj. 

.mate HC 5 — u s \\ La0 (p\ U B^, by the definition of ( s there holds for a; G f2\U-B*. 



y p j i-u s rS) 



|x — SEj I 
-SnS\ 



^ 1 + S a R 5 i v 



and since \x — £i| > 5 a iq in Q, \ UB^, arguing as in estimate (|3.3| gives 



< ^±^Ml + C S 2a - 3 sup R 5 Jl+u d supR«) 

< C(6 a + 5 2a - 3 sup R s A(l + u 5 sup R*), 
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thus, 



i«(n\uB«) 



< C(S a + S 2a - 3 supRt + S 2a - 3 - £ supi?j)(l + u s supR^ 



and finally 

\\u s -u s 

Using the exponent 7, we get 

ll« 5 -fi'lUco ( n\uBf) <C<r(l + 2Bupi2f)(l + « S Bupi2f). □ 

The following lemma gives an estimate for the growth rate of particles in accor- 
dance with the reaction-controlled Lifshitz-Slyozov- Wagner theory. 

Lemma 3.4. For any time t <E (0, T) and 7 > 0, for the growth rates of particles 
holds 



1 

r! 



< C<F(1 + u s supR.f)(l + (1 + (5 7 sup i?f)(l + 2 sup i?f)). 



Proof. Let wf be the capacity potential of the ball Bf with respect to a larger ball 
Bf := B{x h X8 a Rf) for A > 1, i.e., let w t solve 



-Awf = in B^\Bf, ' 
wf = on dB£ 5 , 
wf = 1 in Bf. 



An explicit calculation gives 



wf= 1 



1 - A 



1 - 



XS a Rf 



\X Xi 



and also 



Vwf • n = 



Vtu* • n = 47r — — - 
gB xs 1 — A 



■a d<5 



S a R 



(3.5) 



(3.6) 



(3.7) 



Using equations (|2.ip . (|2.2p . (|2.3p , and the Neumann boundary condition, along 
with the above properties of wf, and integrating by parts, gives 



4n6 2a (Rf) 2 Rf = f V 

JdBf 



u s ■ n 



wfVu 3 ■ n 



IdBf 



VwfVu 6 



u s Vwf ■ n 



dBf 



IdBf 



u s Vwf ■ n 



dB> 



>dB> 



■ 4-7T- 



5 a R 5 



1-A 



R'" 



Rf - u s 



(u s - u s )Vwf 
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where in the last equation we used (|3 . T[) and added and subtracted u s . After 
rearranging, we have 



m - u° - — F 



?A<5 



A 1 1 4Tr\S a RfJ aBt 



A 11 11 47rA(5 Q i? 



,713 



*5 



< + \\u d - u a || io . (nVuB ,), (3.8) 

where in the last step we again used equation (|3.7p . But by using equation (|2.3p 
for on <9i?f we have 

R*\R*\ < 1 + R* \u s \ < 1 + Rf(\\u s - S 5 || Loo( , AuBf) + (3.9) 

Substituting back in (|3.8[) and using Corollary 13. 31 gives the final estimate. □ 

The next lemma ensures that the bounds in the approximation and the growth- 
rate estimates arc indeed uniform. 

Lemma 3.5. For any time t € (0,T), the mean field and the radii of the particles 
are uniformly bounded, i.e., 

u s (t) < C and sup .Rf(t) < C. 



Proof. For the mean field (|3.ip holds 

*' = E / ? TW " sup * (1 + §a sup * f } ? ^/ )E 

and since by Holder's inequality 

i ' i ^ i ' i ^ ^ ' i 

conservation of the total volume of particles gives 

u s < CsupRf(l + 6 a supflf). 

or, using the exponent 7, 

u s < Csupi?f(l + <Fsupijf). (3.10) 

Consider now the set 

A := It I supflf(i) < 



<$t/ 4 J ' 

then, for times t € A, plugging (|3.10[) in estimate (|3.9[) and using Corollary 
gives 

±(Rf) 2 <Csup(i?f) 2 + C. 

Integrating over the time interval (0, T), Gronwall's inequality implies that 

sup, sup teAn[0 T] (i?f) 2 < C(T), 

therefore, [0, T] C A, i.e., the radii are bounded up to time T as is the mean field 
by estimate (|3"TTu) . □ 

Finally, the following lemma gives control over the growth rates of vanishing 
particles and will prove useful for some regularity considerations in the next section. 
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Lemma 3.6. For any time t £ (0, T) such that Rf(t) < l/4sup t $u s (t) and for 
sufficiently small 6, there holds 

— % < R 5 i < — ^"T<0 



Proof. For 



Vt~t < Rf < 2Vt~ 

5 a R s „ S a 



' l + S a Rf\x-x t y 

it can be verified that the function u s — g satisfies the assumptions of the maximum 
principle in Lemma T3. II for the constants a = l/47T(5 2a (i?f ) 2 , thus yielding u s > g 
in ft \ UBf. But since u 5 = g on the boundary dBf, monotonicity implies that 
\7u 5 ■ n > Vg ■ n on dBf and taking the average integrals over dBf we have 

f f s > 1 > 2 

l ~ Rf(l + 6 a Rf)~ Rf 



Moreover, Lemma [3T41 gives 

R-i < u S - -\ + C<F(1 + u s sup Rf ) (1 + (1 + <5 7 sup Rf)(l + 2 sup J?? )) . 

Using now the assumption that Rf < 1/4 sup t ,5 u s and since from Lemma 13.51 it 
follows that for sufficiently small S the 0(<5 7 ) term is uniformly bounded by l/4Rf, 
we get 

1 - 4Rf Rf ARf ~ 2Rf 

Let now 



yi ■= vU — t, V2 ■= 2vU - t 

be sub- and supersolutions that respectively solve 

1 2 
V\ = -7T- . 2/2 = ■ 

2yi y2 



By comparison, we get the lemma's second assertion, i.e., y\ < Rf < y%. □ 

4. HOMOGENIZATION 

In order to pass to the homogenization limit of infinitely-many particles, we need 
first describe the particle-radius density in the limit. To that end, define at any 
time t £ (0, T) the empirical measure vf as 

(0, v\) = J 0(t, Rf(t)) Avf := 1- 4>{t t Rf(t)) for e C c , 

i.e., for functions (f>(t, R) continuous and compactly supported in the radius variable. 
Using now the estimates from the previous section, we can prove the following 

Lemma 4.1. For a subsequence 5 — > and for a function u £ W l,p {Q,T), for 
p < 2, holds 

u 5 ->u mL 2 (0,T), 



u 5 



L 2 {0,T;H\n)). 
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Furthermore, the measures v\ converge to a family v t of 'probability measures such 
that 

J </)dvf — > J (f>a(t)&vt uniformly in t, 
where a(t) denotes the percentage of active particles in the limit. 
Proof. As a consequence of Lemma 13.61 we have 

sup\\Rf\\ LP ^ T) <C(p) forp<2, 
thus, conservation of volume and boundedness of the radii give for p < 2, 

< Csupll-RfH^o.T) < C. 

Lv(0,T) 

Therefore, u s £ W 1,p (0,T), for p < 2, and the compactness following from the 
Rcllich-Kondrachov theorem gives that u s converges to a limit u in L? . 

Taking into consideration that u s is extended to the whole of SI and using the 
lemmas in the previous section, C* 5 converges to u in L 2 (Q) and u s — C" 5 converges 
uniformly to as 5 — > 0, therefore, u s converges to u in L 2 (f2). By the energy 
equality in Lemma 12.11 we have further control over ||Vu <5 ||i2 and thus, we have 
strong convergence in L 2 (0, T; 

For the measures vl holds 

H^ll -supn^^K^^l^l 

in the norm of (C c )*, so for a subsequence 5 — ► there holds that i/f converges 
weakly-* to u t . Furthermore, for positive functions <fi the limit measure v t is non- 
negative and from this it follows that v t becomes zero if there are no particles left 
in the system. 

Choosing now a function tp(t) that depends only on time, we calculate 

The ratio N/Ni is the percentage of active particles at time t. This ratio is bounded 
by 1 and decreasing, therefore it is uniformly bounded in the space BV(fi, T) and 
by the compact embedding of BV(0, T) n L°°(0, T) in L 2 (0, T), it converges in L 2 , 
for a subsequence 6 — > 0, to a limit a £ BV(0, T). 

If we project now the measure v t to the interval [0, T], we get that the projection 
satisfies 

P ro j[o,T] vt = a(t) dt 

and according to Ch. 1, Thm. 10], the decomposition and convergence to v t 
follow from the slicing of measures. □ 

We conclude with the following theorem which states that the limit measure v t 
satisfies the Lifshitz-Slyozov-Wagner equation in a weak sense. Note that in the 
theorem's statement, the initial condition is defined as 

J d>{t,Rf(t)) d*4 := ^r5>(0,itf(0)). 

Theorem 4.2. The measure v t satisfies the Lifshitz-Slyozov-Wagner equation in 
the sense that 

Ji^(t,R)+(u-^j^(t,R)^a(t)du t + y^(0,i?)d^ = 0, (4.1) 
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for all smooth and compactly supported functions <j) € C£°([0, T) x R+), where the 
mean field u is given by 

r Rdv t I j R 2 dv t . 

Proof. We begin by computing the mean-field limit u(t). For a continuous function 
<j>(t) there holds, by the definition of u s , 

?2 1 p2 



J v l + S a R 1 Ni ^ 1 + S a 



1 



Ri 

Ri 



+ 5 a R t 

R 



■dv°. 



l + 5 a R 

Taking the limit S — > on both sides. Lemma |4~T1 gives 

u = J Rdv t J R 2 dv t . 

Consider now a smooth and compactly supported function <p as in the theorem's 
statement. Then, the fundamental theorem of calculus and Lemma 13.41 give 

o = iy 4>(t,Rf(t)) dv & t + J 0(0,^(0)) dv 5 

^4 > (t,Rl{t))+Rl(t)^(t,R!(t))\d^ + J 0(O,i2f(O)) du s 
| ^(t, Rf(t)) + (u s j^j ^(t, Rl(t)) | 6x4 + 0(<F) 
0(O,i?f(O)) d^ 5 . 

The result follows by taking the limit for a subsequence S — ► and using the strong 
convergence of m 5 . □ 

As a concluding remark, we note that the well-poscdncss (existence, uniqueness, 
and continuous dependence on initial data) of the weak formulation (|4.1[) can be 
treated by the methods developed by Niethammer and Pcgo in [14] and [15] . 
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